Holographic Mean-Field Theory for Baryon Many-Body Systems 
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We propose a mean-field approach to analyze many-body systems of fermions in the gauge/gravity 
duality. We introduce a non-vanishing classical fermionic field in the gravity dual, which we call the 
holographic mean field for fermions. The holographic mean field takes account of the many-body 
dynamics of the fermions in the bulk. The regularity condition of the holographic mean field fixes the 
relationship between the chemical potential and the density unambiguously. Our approach provides 
a new framework of gauge/gravity duality for finite-density systems of baryons in the confinement 
phase. 
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Many-body physics of strongly-correlated fermions is 
one of the central subjects in modern physics. Strongly- 
correlated fermions appear in various places, such as the 
systems of cold atoms, high-T c superconductors, neutron 
stars and the quark-gluon plasma. However, the the- 
oretical description of many-body systems of strongly- 
correlated fermions is still a challenge owing to their 
non-perturbative natures. One of the successful compu- 
tational frameworks for strongly-interacting fermions is 
the lattice gauge theory, but the Monte Carlo simulation 
at finite baryon density suffers from the sign problem. 

Recently, a great deal of attention has been paid to 
the AdS/CFT correspondence [iHll because of its abil- 
ity for non-perturbative computations and the absence 
of the sign problem. However, the current framework 
of the AdS/CFT correspondence for finite baryon den- 
sity still has room for improvement in the confinement 
phase. It can be summarized as follows: I) dynami- 
cal baryons at finite density have not been satisfacto- 
rily incorporated, and 2) the Gubser-Klebanov-Polyakov- 
Witten (GKP-Witten) prescription @, 0| to relate the 
baryon chemical potential to the baryon density has an 
ambiguity Q. 

In this paper, we shall propose a new approach, a holo- 
graphic mean-field theory for fermions, to resolve these 
problems. The self-consistent incorporation of dynamical 
baryons in our approach provides an eigenvalue equation 
which fixes the relationship between the baryon chemical 
potential and the baryon density unambiguously. 

Let us revisit the above mentioned problems in detail 
below. The quarks are realized as fundamental strings 
attached to the flavor brane in the gravity dual [H[ . Each 
end point of the strings carries a unit charge with re- 
spect to a U(l) gauge field on the flavor brane, and 
the quark density can be evaluated from the "electric 
field" in the radial direction at the boundary through 
the Gauss-law constraint. The quark chemical potential 
is thus the time component of the U{1) vector poten- 
tial, Aq, at the boundary [6], which is conjugate to the 
electric field. This fits the conventional GKP-Witten pre- 
scription: the boundary value of Aq gives the source (the 



chemical potential); the expectation value of the conju- 
gate operator (the density) is given by the derivative of 
Aq at the boundary with respect to the radial direction. 
These quantities are boundary conditions for the Gauss- 
law constraint which may be chosen arbitrarily. However, 
another physical condition in the bulk, i.e., the bulk con- 
dition, fixes the ambiguity in the relationship between 
them whenever the GKP-Witten prescription works. To 
the best of the authors' knowledge, only the example 
where the bulk condition is yet to be known is the finite 
density systems of baryons in the confinement phase. 

The origin of the ambiguity can be understood in the 
following way. The global U(1)b (Z7(l)-baryon) symme- 
try is enhanced to the U(l) gauge symmetry in the grav- 
ity dual. The constant shift of the boundary value of Aq 
can be a gauge transformation which does not alter the 
physics in the gravity dual whereas the constant shift 
of the chemical potential has definite physical meaning 
in the boundary theory. This contradiction is resolved 
when the flavor brane intersects the horizon: Aq has to 
vanish at the horizon to make the one-form well-defined 
there @, Q ■ This removes the ambiguity of the constant 
shift of the boundary value of Aq. However, we cannot 
employ this prescription if the flavor brane does not in- 
tersect the horizon: we need another bulk condition from 
which the ambiguity in the constant shift of Aq is elimi- 
nated. 

More serious problem is the difficulty to incorporate 
many-body systems of dynamical baryons. In the phase 
where the flavor brane is away from the horizon, quarks 
are confined. The U(\)b charge is carried by baryons in 
this case. We assume that the number of the colors, N c , 
is odd: the baryons are fermions. A single baryon is holo- 
graphically realized as a baryon vertex [9( which consists 
of N c fundamental strings and a D-brane wrapped on a 
compact subspace. This is a fermionic object that lives 
inside the bulk of the dual geometry. 

In order to prepare a system of finite density of 
baryons, we need to deal with the many-body physics 
of the baryon vertices. We expect that their distribu- 
tion along the radial direction in the dual geometry is 
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non-trivially determined by their dynamics. However, it 
is very difficult to solve the many-body dynamics of the 
baryon vertices within the current technology of super- 
string theory: we need a good approximation. Although 
they can be described as solitons (or Skyrmion-like ob- 
jects) [lol - H^ I on the flavor brane, the many-body physics 
of solitons on the curved spacetime is still hard to solve 
[l3| : we need further simplification. In most of the lit- 
erature 14[ , many-body systems of baryons in the grav- 
ity dual have been constructed just by superposing the 
baryon vertices at the same location in the radial direc- 
tion. However, the dynamics of the baryon vertices is 
not fully taken into account in this method: the system 
is over-simplified. 

One way to overcome the difficulty is to approxi- 
mate the baryon vertex as a point-particle fermion rather 
than the wrapped D-brane with N c strings, and to de- 
scribe their many-body dynamics in terms of a Dirac 

field m h m on the dual geometry. Indeed, this ap- 
proximation makes sense at the stron g co upling where 
the size of the soliton becomes small [111, [l2j. In this 
paper, we proceed along this approximation. We con- 
struct a system of finite baryon density by using the non- 
vanishing bulk Dirac field, which we call the holographic 
mean field for baryons. 

Since the baryon vertex carries the U(l) charge, the 
holographic mean field is also requested to carry it. Then 
the mean field induces non-zero expectation values of the 
fermion bi-linear operators, such as the charge density 
which is detected by the derivative of the U(l) gauge 
field at the boundary. The dynamical distribution of the 
baryon vertices along the radial direction is represented 
by the configuration of the mean field. Furthermore, the 
coupled equations of motion for the mean field and the 
U(l) gauge field give an eigenvalue equation from which 
the chemical potential is unambiguously fixed. It will also 
be shown that the chemical potential at the zero-density 
limit automatically agrees with the mass of the baryon. 
Thus, the problems 1) and 2) are resolved in harmony in 
our approach. 

Let us consider a five-dimensional (5D) dual geometry 
without an event horizon whose metric depends only on 
the radial direction, for example. It is convenient, al- 
though not necessary, to make a coordinate transforma- 
tion after which the geometry is given by the conformally 
flat metric, 
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H{w)(ji^ v dx' l dx v - dw 2 ) , 



(1) 



where 77 MI/ — diag(+, —,—,—) and w is the radial direc- 
tion. An example of H(w) is given in 11] . We introduce 



a 5D Dirac spinor ^> which interacts with a bulk £7(1) 
gauge field Am on this geometry. The actions for them 



5* = J d A xdw [Wr M (d M - iqA M )V - m 5 (ttf)#*] , 
Sa = I d i xdwCA , (2) 



where q and 772,5(10) denote the charge and the 5D mass 
of the fermion, respectively. T M are written in terms 
of the four dimensional (4D) Dirac matrices as T v = 7^ 
= 0,1,2,3) and iT w = 7 5 . The possible boundary 
terms and the counter terms are omitted here. Note that 
the Dirac field is on the curved spacetime ([TJ, but the 
vielbein and the spin connection are absorbed by a suit- 
able field redefinition of the fermion field [HI]. The w 
dependence of 7715(10) reflects the curvature of the geom- 
etry. Typically, 771.5(10) diver ges at the boundary, and its 

We assume that Sa is 
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explicit example is given in 
a functional of On Am, which may be given by a Dirac- 
Born-Infeld (DBI) action on the flavor brane, for exam- 
ple. However, its explicit form is not important at this 
stage. We employ the probe approximation such that 
the fermion field and the gauge field do not modify the 
geometry: we regard S = S^ + Sa as the total action up 
to the constant action of the gravity sector. 

We now explain how we formulate the holographic 
mean-field approach. We work in the A w {x,w) = 
gauge, and replace the gauge field with the mean field 
as 

Aq(x,w) ->• A (w), Ai(x,w) -> Ai(w) , (3) 
where i = (1, 2, 3). We impose 

(4) 



■^■0 (^0 I boundary ^ ' 



Ai(w)\, , =0. 

v ' I boundary 



where \i is the chemical potential associated with the 
fermion charge. We also introduce the mean field for the 
fermion as 



(5) 



The non-vanishing mean field realizes the dynamical dis- 
tribution of the charge in w direction, which is given by 



p(w) = q^(w)ty(w) . 



(6) 



The integration of p(w) over w gives the 4D charge den- 
sity n. In the present analysis, we switch off the fermionic 
source: 



*(w) 



boundary 



= 0. 



(7) 



We also require that the mean fields are regular. 

Now A' = 8 w Aq is given by the Gauss-law constraint, 

d w [d£ A /dA' (w)} - q^(w)*(w) =0, (8) 

and A\ = d w Ai is given by 

d^dCA/dA'^w)} - q^(w)r°r^{w) = . (9) 
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The scalar potential Aq(w) is given by using A' Q (w): 
A {w)= f i + dw'A' {w')=ii + A {w). (10) 

J boundary 

However, the boundary value /i, that is the constant of 
integration, is arbitrary at this stage. In other words, 
whatever value of /x can give the same charge distribution. 
This has been a problem in holographic models at finite 
density. This problem will be resolved by solving the 
equation of motion for ty(w): 

\iT w d w + qT a A (w) + qVA^w) - m 5 (w)] 9(w) = 

(11) 

that depends on the boundary value of Aq (w) explicitly. 

Let us see how it works. We write ^(w) in terms of the 
two-component spinors 1 $>± as *ff(w) = ( V I'+,4'_) T with 
the Dirac matrices given by 



On the other hand, at \i = 0, all the mean fields vanish, 
then (fT4l) is reduced to be 



1 

1 



1 

-1 



a 1 
-a 1 



(12) 

We can take A\(w) = A2(w) = 0, and hence a 3 ^>± = 
*Sf± without loss of generality by virtue of the rotational 
invariance. From (|11[) , we obtain 

[d w - qA 3 (w)]^+ - [m 5 (w) + qA {w)}^- = t?/x*_ , 

[d w + qAsiw)]®- - [m 5 (w) - qA (w)]^ + = -qfi^+ . 

(13) 

Solving (fl~3|) coupled with ([8} and ©, q/j, is obtained as 
the eigenvalue for a given value of n = J dw q^ (w)^ (w) . 
The point here is that the eigenvalue equations (fl3|) are 
nonlinear equations after rewriting the gauge fields in 
terms of ty + and and the amplitudes of ^+ and ^_ 
are not arbitrary if we specify qp,. If we start with a 
given qp, the density n is determined as a function of 
qfi. In this way, the holographic mean-field theory pro- 
vides the equation of state (EOS): the relation between 
the 4D charge density n and the chemical potential a is 
determined dynamically. The new bulk condition to fix 
the ambiguity of /i in the GKP-Witten prescription is the 
regularity of the fermionic mean field under ([7]) . 

Once the mean fields are obtained, the energy spec- 
trum of the physical excitations can be computed by in- 
troducing the fluctuations ij)(x, w),a v (x, w) on top of the 
mean fields as ^(x,w) — ty(w) +ip(x,w) and A u (x,w) — 
A v (w) + a v (x,w), and solving the coupled equations of 
motion for the fluctuations. One of the equations of mo- 
tion is given by 

\iT w d w + r°( Po + qA Q (w)) + T 3 qA 3 (w) + f-p 
-m 5 (w)]ij}(po,p, w) +r l 'ty(w)qa l/ (p ,p,w) =0. (14) 

Let us discuss the zero density limit, n — > +0. (jHJ) and 
||SJ) tell us that A (w) -)■ 0, A 3 (w) -4 0, and (HU) becomes 



d w ip+ 



■ m 5 (w)ip+ 



-mfip+ , 



(16) 



where ip — (ip + ,ip-) T and we used T v p u xj) — mfij}. 

The eigenvalue my of (|16[) for the non- vanishing eigen- 
function under the boundary condition ib\. . = is 

17 ' I boundary 

the mass of the fermion. Therefore, the threshold of the 
chemical potential for the non- vanishing density given by 
(|15p coincides with the fermion mass, i.e., qu — > rrif at 
the limit n — > +0. We do not find any nontrivial solu- 
tions for ty(w) for \qu\ < rrif, where there is no formation 
of the Fermi surface. 

Now we demonstrate how to obtain the EOS for baryon 
many-body systems. We apply our method to the model 
given in [llfl where the baryon field is introduced into the 
Sakai-Sugimoto model 17J. The boundaries in this model 
are located at w = ±w max , thus the boundary conditions 
for the mean fields are Ao(±u; max ) = /J,, A 3 (±w max ) = 0, 
$(±!« ma i) = 0. Here /x is the baryon chemical potential 
and q = 1 . We assume that the theory is invariant under 
the 5D parity (t,x,w) — > (t, — x, — w) [17j. Therefore, 
Aq(w) (A 3 (w)) is an even (odd) function of w, and the 
regularity condition leads to ^4o(0) = (^3(0) = 0). 
From the parity invariance, ty + ( v I / _) is either an even 
(odd) function of w or an odd (even) function. Since 



(Tl3l) is invariant under q) -q), the 

assignment of the parity is a convention up to the sign of 
q. In this demonstration, we assign even parity to 
and the regularity condition we employ is ^+(0) = 0. 
We use the DBI action of the D8-branes as Sa ■ 

S A = -C I d 4 xdwU 1/4 (w) 



(15) 



x ^H*(w){H*(w) {2^ s /N c f[{A' Q {w)Y - (A' 3 (w))*}} , 

(17) 

where C = (27r)- 7 (4ir 2 /3)ls 11/2 X 1/i (2M KK )- 1 ^N f N c , 
and the A v {w) here is N c times A u (w) in the standard 
convention. Here, Mkk is the energy scale of the theory, 
A is the 't Hooft coupling, Nf is the number of flavors, 
and l s is the string length which does not show up in the 
end of the calculation. The definitions of U(w) and H(w) 
are given in (ll| . 

We fix Mkk and A in the same way as in !17]: Mkk = 
949 MeV and A7V C = 50. We consider N f = 2 case, and 
use the normal nuclear matter density no = 0.16 fm~ 3 to 
scale the baryon number density n. Figure[T]shows the re- 
sultant EOS compared with the one for the free baryons. 
The chemical potential \x increases as the density grows. 
The increase of the chemical potential of free baryons is 
very tiny on the figure, whereas that of our dynamical 
baryons is considerably larger: this can be understood as 
an effect of repulsive interactions among the baryons. 
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FIG. 1: The EOS of the model (solid) compared with the 
one for the free baryons (dashed) given by u. = \/m\ + pp, 
where m n and pf are the mass and the Fermi momentum of 
the baryons, respectively. The horizontal axis is n/no, while 
the vertical axis is p/m n . 
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FIG. 2: The baryon charge distribution p(w)/n for n/no = 
0.1 (solid), 1 (dashed) and 2 (dotted). The boundary is lo- 
cated at w = iu max — 3.64/Mkk- The unit of w on the figure 
is 1/Mkk. 



Figure [2] shows the density dependence of the baryon 
charge distribution p(w), which is indeed regular. The 
distribution shifts to the boundary side as the density 
grows. This implies that taking account of the distri- 
bution of the baryon charge along the radial direction 
becomes more important in the higher density region. 

We expect that our approach is applicable to holo- 
graphic descriptions of various systems of fermions [l8| . 
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Note added: After the first version appeared at the 
arXiv, we were informed on a recent related work [l9| 
where the Fermi distribution of the boundary theory is 
introduced as an input to obtain the EOS. 
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